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Abstract 

This paper describes an algorithm for selecting a consistent set within 
the consistent histories approach to quantum mechanics and investigates its 
properties. The algorithm select from among the consistent sets formed by 
projections defined by the Schmidt decomposition by making projections at 
the earliest possible time. The algorithm unconditionally predicts the possible 
events in closed quantum systems and ascribes probabilities to these events. A 
simple random Hamiltonian model is described and the results of applying the 
algorithm to this model using computer programs are discussed and compared 
with approximate analytic calculations. 

1 Introduction 

It is hard to find an entirely satisfactory interpretation of the quantum theory of 
closed systems, since quantum theory does not distinguish physically interesting 
time-ordered sequences of operators. In this paper, we consider one particular line 
of attack on this problem: the attempt to select consistent sets by using the Schmidt 
decomposition together with criteria intrinsic to the consistent histories formalism. 
For a discussion of why we believe consistent histories to be incomplete without a set 
selection algorithm see |l],[2|] and for other ideas for set selection algorithms see . 
This issue is controversial: others believe that the consistent histories approach is 
complete in itself [0-0]. 
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1.1 Consistent histories formalism 



We use a version of the consistent histories formalism in which the initial conditions 
are defined by a pure state, the histories are branch-dependent and consistency is 
defined by Gell-Mann and Hartle's medium consistency criterion eq. (0). We restrict 
ourselves to closed quantum systems with a Hilbert space in which we fix a split 
Ti. = 7ii <8> 71-2', we write dim(TCj) = dj and we suppose that d\ < d 2 < oo. The 
model described in sec. |2] has a natural choice for the split. Other possibilities are 
discussed in ||. 

Let be the initial state of a quantum system. A branch- dependent set of 
histories is a set of products of projection operators indexed by the variables a = 
{a n , a n _i, . . . , q;i} and corresponding time coordinates {t n , . . . , ti}, where the ranges 
of the ak and the projections they define depend on the values of ak-i, ■ ■ ■ , on, and 
the histories take the form: 

C a = f£(f n ; . . . , ai)i*£i(^ i; a n _ 2 , ...,«!)... P^(h) . (1) 

Here, for fixed values of aik-i, • • • , ai, the P% Atk] otk-i, • • • , «i) define a projective 
decomposition of the identity indexed by a^, so that X) a PaA^k] a k-i, ■ ■ ■ , a i) — 1 
and 

Here and later, though we use the compact notation a to refer to a history, we intend 
the individual projection operators and their associated times to define the history. 
We use the consistency criterion^ 

D aP = 0, Va ^ 0, (3) 

which Gell-Mann and Hartle call medium consistency, where D a p is the decoherence 
matrix 

D af3 = Tr(C aP Cl). (4) 

1 For a discussion of other consistency criteria see, for example, refs. |Io|-|l3[] . 



Probabilities for consistent histories are denned by the formula 



p(a) = D aa . (5) 

With respect to the H = Hi <8> H2 splitting of the Hilbert space, the Schmidt 
decomposition of \ip(t)) is an expression of the form 

|V(i)) = E[Pi(t)] 1/2 k(*))i®K(*)}2, (6) 
1=1 

where the Schmidt states {\wi)i} and {1^)2} form, respectively, an orthonormal 
basis of Hi and part of an orthonormal basis of H2, the functions Pi(t) are real and 
positive, and we take the positive square root. For fixed time t, any decomposition 
of the form eq. (|(J) then has the same list of probability weights {pi(t)}, and the 
decomposition ([]) is unique if these weights are all different. These probability 
weights are the eigenvalues of the reduced density matrix. 

The idea motivating this paper is that the combination of the ideas of the consis- 
tent histories formalism and the Schmidt decomposition might allow us to define a 
mathematically precise and physically interesting description of the quantum theory 
of a closed system. We consider constructing histories from the projection operators^ 

P i (t) = \w i (t)) 1 (w i (t)\ 1 ®I 2 and P = h®I 2 -EiPi(t), ( 7 ) 

which we refer to as Schmidt projections. If dim7ii = dimH2 the complementary 
projection P is zero. In developing the ideas of this paper, we were influenced 



in particular by Albrecht's investigations |L4|,|l5j of the behaviour of the Schmidt 
decomposition in random Hamiltonian interaction models and the description of 
these models by consistent histories. 

2 A Random Hamiltonian Model 

Consider a simple quantum system consisting of a finite Hilbert space H = Hi ® H2 
(dim?-^j = dj), a pure initial state \ip(0)) and a Hamiltonian drawn from the GUE 



2 Thcre are other ways of constructing projections from the Schmidt decomposition though 
for the model considered in this paper the choices are equivalent. 



(Gaussian Unitary Ensemble), which is defined by 



P(H) = N- 1 exp{-Tr[(A# + /i) 2 ]}, (8) 

where N is a normalisation constant. 

The GUE is the unique ensemble of Hermitian matrices invariant under U(d) 
with independently distributed matrix elements, where d = dimH = d\d 2 . The 
GUE is also the unique ensemble with maximum entropy, — / dH P(H) log P(H), 
subject to E[Tr(H)] = fj, and E[Tr(H 2 )] = A. The book by Mehta JIJ] contains a 
short proof of this as well as further analysis of the GUE and related ensembles. All 
the results concerning the GUE in this thesis can be found in this book or in the 
appendix. 

This model is not meant to represent any particular physical system, though 
Hamiltonians of this from are used in models of nuclear structure and have often 



been studied in their own right (see [16 , [T7j and references therein), and a large class 



of other ensembles approximate the GUE in the large d limit. 

Because H is drawn from a distribution invariant under U(d) there is no preferred 
basis, no distinction between system and environment degrees of freedom and no 
time asymmetry. In other words the model is chosen so that there is no obvious 
consistent set: we do not already know what the answer should be. Moreover it 
does not single out a pointer basis that one might associate with classical states, 
so that the Copenhagen interpretation cannot make any predictions about a model 
like this in the t —>■ oo limit. If an algorithm works for this model, when there are 
no special symmetries, it should work for a wide variety of models. The question 
whether a pointer basis can arise dynamically using Schmidt states was addressed 
by Albrecht in |]1^ , [15| , but no general prescription emerged from his study. Albrecht 
also studied the relationship between Schmidt states and consistent histories, and 
his studies suggested that the relationship was complicated. 

The model considered here generalises Albrecht's model: the Hamiltonian for 
the entire Hilbert space is chosen from a random ensemble. Albrecht also used a 
different distribution, but as we explain below the GUE seems more natural, though 
this it probably makes little difference. 

Without loss of generality, we take \i = and A = 1/2. With this choice and 



using the Hermiticity property of H eq. (^) becomes 



P(H) = n 



-n 2 /22-n/2jy-l 



Hu\ 2 /2 



(9) 



i<j i 



Therefore the diagonal elements are independently distributed, real normal random 
variables with mean and variance 1 and the off-diagonal elements are independently 
distributed, complex normal random variables with mean and variance 1. 

Since the Hamiltonian is invariant under U(d) the only significant degrees of 
freedom in the choice of initial state are the initial Schmidt eigenvalues (the eigen- 
values of the initial reduced density matrix.) The usual choice in an experimental 
situation is an initial state of the form = \u)i ® \v)2 which corresponds to a pure 
initial density matrix. A more general choice in the spirit of the model is to draw 
the initial state from the U(d) invariant distribution subject to fixed rank n. This 
is equivalent to choosing the first n eigenvalues to be components of a random unit 
vector in R n and the remaining d\ — n components to be zero. 



The calculations in this section are an attempt to gain insight into the expected 
properties of prediction algorithms applied to the random model. These calcula- 
tions rest on a large number of assumptions and are at best approximations, but the 
conclusions are borne out by numerical simulations and the calculations do provide 
a rough feel for the results that different algorithms can be expected to produce. In 
particular they suggest that there are only narrow ranges of values for the approx- 
imate consistency parameter which are likely to produce physically plausible sets 
of histories. These calculations may also be applicable to other models since this 
model makes so few assumptions and the interaction is completely general. 

In a random model there is no reason to expect exactly consistent sets of histo- 
ries formed from Schmidt projections to exist, so only parameterised approximate 
consistency criteria such as the frequently used criterion plB|.|T9[| 



3 Analysis 



D aP \<e(S), 



(10) 



or the DHC (Dowker-Halliwell Criterion) 



\Re(D a p)\ < e(D aa D^y/ 2 , Wa^P, (11) 

or 

\D a/3 \ < e{D aa D pp ) l '\ Va=£p. (12) 

are considered in this paper and e will always be the consistency parameter in these 
equations. We shall only discuss medium consistency criteria: the results for weak 
consistency are qualitatively the same. 

Approximate consistency criteria were analysed further in ref. As refs. 

fT2j . p~3|j explain, the DHC has natural physical properties and is well adapted for 
mathematical analyses of consistency. We adopt it here, and refer to the largest 
term, 

max{ iDapKDaaDpp)' 1 / 2 : a, f3 G S , a ^ (3 , and D aa , % ^ } , (13) 

of a (possibly incomplete) set of histories S as the Dowker-Halliwell parameter, or 
DHP. 

If an absolute approximate consistency criterion is being used there are strong 
theoretical reasons for imposing a parameterised non-triviality criterion JTB|. How- 
ever, if the approximate DHC is being used one is not needed, though it is convenient 
to introduce one for computational reasons. The non-triviality parameter (which we 
shall always write as 5 in this chapter) can be taken very small if the DHC is used 
and is not expected to influence the results — except possibly for the first projec- 
tions — and the numerical simulations show that this is indeed the case. We shall 
refer to histories with probability less than or equal to 5 (relative or absolute) as 
trivial histories and a projection that gives rise to a trivial history as a trivial pro- 
jection. There are no absolute reasons for rejecting set of histories containing trivial 
histories — if 5 is sufficiently small and there are not too many they are physically 
irrelevant — though obviously sets are preferable if all the histories are non-trivial. 
However, an algorithm must produce results that are approximately the same for 
a range of parameter values if it is to make useful predictions, and trivial histories 



will almost certainly vary as 5 is changed. If the DHC is used, generically all the 
later projections will also change, since trivial histories can significantly influence the 
consistency of later projections. If an absolute consistency criterion is used trivial 
projections are more likely to be consistent than non-trivial projections so for many 
values of the parameters only trivial will be projections are made. 



3.1 Repeated projections and relative consistency 

Consider a history a extended by the projective decomposition {P, P} and the fur- 
ther extension of history P\a) by {P(t), P{t)}. This was discussed for the DHC in 



refs. [y,[13[ and the DHP for this case was shown to be 

\(a\PP\a)\ 



(14) 



||P|a)|| ||PP|a)|| 

The reprojection will occur unless e, the approximate consistency parameter, is 



smaller than (|i4|). It is easy to show that the time evolution of Heisenberg picture 
Schmidt projections is 

P = i[H-B®I,P], (15) 

where H is the Hamiltonian, 

B=Y. QkPr _ Qm , (16) 

kj L m Prn — Pk 

Qk are projection operators (in H{) on to the Schmidt eigenspaces, pk their respective 
(distinct) eigenvalues and p r the derivative of the reduced density matrix. 

In analysing ( |T4"D and similar expressions we make the following assumptions. 
First that \a) is uncorrelated with the Schmidt states — this generally is a good 
approximation when there are a large number of histories. Second that B ® / is an 
operator drawn from the GUE with unspecified variance independent of the other 
variables — in some situations this assumption is exact but it generically is not. 

Let G = H — B ® I , an element of the GUE with variance a, then using eq. ([15]) 
(Ml is 

i/™ i pnp\n,\\ 

(17) 
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Because G is drawn from a distribution invariant under U(d) and is independent of 



P\ot) and P\a), (|17|) can be simplified by choosing a basis in which P|a)/||P|a)|| = 
(1,...,0) and P|a)/||P»|| = (0,1,..., 0). (0) becomes 









[J2r>k>l 


Zk 


2] 1/2 



where r = rank(P) and Zk = Ca(fc+i). Since {Gij,i < j} is a set of independent, 
complex, normal random variables, flT8|) is the square root of a B(l, r — 1) random 
variable^]. 

Suppose we choose e so that reprojections will occur with some small probability 
q — note that only choosing e = will definitely prevent all repeated projections. 
The probability of (|i8|) being less than e is 

l-(l-e 2 ) r -\ (19) 



Therefore if 



[1 _ (1 _ g )l/(r-D]l/2 



log(l - q) 



1/2 



+ 



log(l - q) 



3/2 



<d-VV-l g(l- g ) 



(20) 



a reprojection will occur with probability ~ q. 

However, it is shown in || that the DHC cannot prevent trivial reprojections on 
the initial state if the initial density matrix has less than full rank. If the initial 
density matrix has rank one then the first projection will always be made with 
probability 5. A non-triviality criterion can then work in conjunction with the 
DHC to prevent further trivial extensions. Suppose either that the initial density 
matrix has rank greater than 1 and P n and P m are two projections onto the non-zero 
eigenspaces, or assume that the rank is one and P n is a projection onto the initial 
state and P m is a projection making a history of probability 5. In either case, let P& 
be projection onto the null space. To prevent the trivial projection P& being made 



3 B(p, q) :- a beta random variable with parameters p and q. This has a density function 
oc t p_1 (l — t) 9 ~ l . A B(l,r-1) random variable has the same distribution as that of the inner- 
product squared between two independent unit vectors in C r . 



the parameters 5 and e must be chosen to satisfy || 



V5\{^\P m P^PM\ > e\\P m m\ ||P fe P n |^>|| 2 . (21) 

Though the probability distribution for this is complicated, the approximate relation 
between S and e can be estimated by squaring both sides of eq. ( 21 ) and then taking 
the expectation. Note that treating G as an element of the GUE is exact in this 
case as the terms involving B ® I are identically zero. Using results from eq. ([3(|) . 
eq. (|2T| ) becomes 

<5> e 2 (r + l)||P n |^>|| 2 , (22) 

where r is the rank of Pk- By assumption ||P n |^)|| is order one and r < d, so if 
5 > de 2 initial reprojections will not occur. The results are the same for a relative 
non-triviality criterion since instead of eq. (^) we have S > e 2 (r + 1). 

3.2 Repeated projections and absolute consistency 

An algorithm using an absolute parameterised consistency criterion will make noth- 
ing but trivial projections unless a parameterised non-triviality criterion is also used, 
so only algorithms with a non-triviality criterion are considered. 

Let t e denote the latest time that the reprojection is approximately consistent 
and ts the earliest time at which the extension is absolutely nontrivial. We see from 



refs. [y,[T3| that, to lowest order in t, 



t 5 = VdWPPla)]]- 1 (23) 
t £ = eKalPPPla}]- 1 . (24) 



ts > t e implies 



VS\(a\PPP\a)\> e\\PP\a)\\. (25) 
Again we choose e so that reprojections occur with probability q and assume that 



|P|ct)|| and ||P|o;)|| are order one, so that eq. (f25|) can be written 

e/y/6. (26) 













2] 1/2 



The l.h.s. is the same random variable as in eq. (|18D so 5 and e must be chosen so 



that 

e/VS = rf^V- log(l - q) ■ (27) 

The assumption that ||P|a)|| and ||P|a)|| are order one will obviously not always be 
valid. As more projections are made the probabilities of the histories will decrease. 
When both probabilities are 5 eq. ( |25|) is 









[J2r>k>l 




2] 1/2 



so t/5 = cP 1 / 2 '\j — log(l — q). If reprojections of smaller probability histories are to 
be prevented this choice of parameters is clearly more appropriate than eq. fl2"7p. 

This analysis has picked a very conservative upper bound for e to prevent re- 
peated projections, since decoherence matrix terms with the other histories will 
tend to reduce the likelihood of repeated projections, and thus allow larger values 
of e to be used. A more detailed analysis suggests that for relative and absolute 
consistency r can be treated as much smaller than d so that choosing 5 a small 
factor larger than e 2 or e respectively, are sufficient conditions. 



3.3 Projections in the long time limit 

The previous subsection has shown how e and 5 affect the probability of repeated 
projections: this subsection calculates how they affect the probability of projections 
as t — > oo. In infinite dimensional systems, off-diagonal terms of the decoherence 



matrix for quasiclassical projections often tend to zero as t increases |[20|| . In the 
limit d — > oo one would also expect this for Schmidt projections in this model - 
though the limit only exists for initial density matrices of finite rank. 

Consider the DHP for a Schmidt projection extending history a from the set 
of normalised exactly consistent histories {\a),\Pi),i = 1, . . . , k} as t — > oo. For 
t — > oo and for large k the Schmidt states are approximately uncorrelated with the 
histories. The DHP for an extension {P, P} of history a is 



™* no N > i4 i ' = ■ (29) 



m\p 
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a)\ 


\\p 


W)\\ 



Since (fii\a) = for all i, eq. ( p9[) is equal to (within a factor of a/2) 

|(ft|P|g)| 

max — — — _ — — . (d(JJ 



The cumulative frequency distribution for (|30D squared is calculated in ||] as 



max 
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a)ll 2 



<*1 =E(- i rf*l( i -"^) ,i " 1 wi, (si) 



which approximately equals [1 — e _<iA ] fc when dk 2 X 2 = o(l). This is the probability 
that a pair of projections acting on one history in a set of k + 1 consistent histories 
satisfies the medium DHC with parameter e = y/\. There are n p = 2 mm ( <il ~ 1 > ci2 ) — 
1 distinct choices for the projections so the DHP (to within a factor of v2) for 
extending a with any Schmidt projection is 



mm max - — ^_ , (32) 
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«>ll 



where {Pj, P^} range over all n p binary partitions of the basis Schmidt projections. 
The distribution for this random variable is hard to calculate but if we assume that 
the DHP's for each {Pj, Pj} are independent the cumulative distribution function is 

1 - {1 - [1 - e - dX+oW ] k } n " . (33) 



This assumption is obviously very approximate since the different projections are all 
formed using the same basis. However, treating the {Pj, Pj} as independent in ( |3"2"| ) 
will be a lower bound for the exact result and ( j30|) will be an upper bound for the 
exact result. 

Suppose now we wish to choose e so that the probability of making a projection 
at a large time is p, where p is close to one. Then from eq. (|33|) 

e 2 = -l/dlog{l - [1 - (1 - p) 1 ^"] 1 ^} . (34) 



For large k 



e 2 w -l/dlog{-l/Mog[l - (1 -p) 1/np }} w l/dlog(fc) 



(35) 



This calculation has involved a lot of assumptions and approximations, but it should 
accurately reflect the behaviour for d > ^ 2 > 1. The logarithmic dependence on 
k is a generic feature of extreme order statistics and so is the independence of the 
answer from the other factors p and n p . The 1/y/d dependence is also expected 
because the mean value of the square of the inner product between two random 
vectors in a d dimensional space is l/d. 

A generic feature of asymptotic extreme order statistics such as the previous 
calculation, is the slow rate of convergence and this calculation is only expected 
to be accurate for very large d and k. For actual application to particular models 
the exact distribution can be calculated by a computer program using Monte-Carlo 
methods. 

For each sample the program randomly picks a set of exactly consistent histories, 
and then calculates the DHP for all n p combinations of projections on one partic- 
ular history. Ten thousand samples were sufficient to produce smooth cumulative 
frequency distributions. These can be inverted to calculate e(p,k) as in fig. (JJ), 
where for example e(.99, k) is the solid curve. 

The same arguments apply for absolute consistency as t — » oo, but since the 
consistency requirement is not normalised the expected DHP values will be reduced 
by a factor of 1/ \fk — the average value for the length of a history when there are 
k histories. 

These calculations suggest that choosing e = 0{l/\fd) is most likely to produce 
histories with a complicated branching structure and many non-trivial projections. 
In the next section we discuss the results of simulations for all values of the param- 
eters and show that they agree with these theoretical calculations. 

4 Computer simulations 

The computer programs are explained and listed in ref. |4]]. The results described 
here were carried out with a system of dimension 3, with an environment of dimen- 
sion 15 and with either medium absolute consistency or medium relative consistency 
(DHC). Fig. ([[]) gives the probability distribution for the DHP plotted via percentile 
curves as a function of the number of histories in the long time limit. For example, 
this graph shows that with e > 0.3 ~ 2/\fd projections will almost certainly be 
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Figure 1: estimating e 



consistent for any number of histories, whereas for e < 0.05 ~ .3/y/d projections 
will probably only be consistent when there are two or three histories. When there 
are twenty histories it shows that for 98% of the time the DHP will be between 
0.15 and 0.25. Fluctuations will only occasionally (1% of the time) lie below the 
solid line and if e has this value (for a particular number of histories) then although 
projections will probably occur they will occur as a result of large fluctuations from 
the mean. Therefore one would expect the projections to occur at widely separated 
times and if e is changed only slightly the times generically to change completely, 
and indeed computer simulations show this. 

The simulations described here were run for ten thousand program steps or until 
thirty histories had been generated. For a given set of parameters, many simulations 
with different Hamiltonians and initial states were carried out and were found gener- 
ically to produce qualitatively the same results, though only individual simulations 
are described here. 

One way to look at the results of a simulation is to look at the probability tree 
associated with the set of histories such as fig. (0a). The root node on the far left 
represent the initial state, the terminal nodes represent the histories and the other 
nodes represent intermediate path-projected states. Each node has a probability and 
the lines linking the nodes have an associated projection operator and projection 
time. The projections associated with lines emanating to the right from the same 
node form a projective decomposition and all occur at the same time. The scaling 
of the axis and relative positions between the nodes is arbitrary, only the topology 
is relevant. For example, in fig. fl^a) the probabilities for the histories are 0.42, 0.25, 
0.05, 0.02, 0.02, 0.01, 0.08 and 0.14 — the probabilities of the terminal nodes. 




(a) probability tree (b) consistency statistics 



Figure 2: example 

Another useful interpretative aid is a graph of the consistency statistics fig. (0b). 
This graph shows the DHP for the most consistent non-trivial extension. At times 
where no Schmidt projections result in non-trivial histories no points are plotted, 
though there are no such times in fig. (0b). The program makes a projection when 
this value is e. The flat line indicates e and the crosses indicate when projections 
have occurred — in this case at times 0, 1, 11, 12, 36 and 65 (approximately). A 
graph of the projection times will also be used sometimes, for example fig. (f|b). 

When any Schmidt eigenvalues are equal their eigenspaces becomes degenerate 
and the corresponding Schmidt projections are not uniquely defined. The reduced 
density matrix varies continuously in this model and it will only be degenerate for 
a set of times of measure zero so generically it is possible to define the Schmidt 
states so that they are continuous functions of t for all t. This was not found to be 
necessary in the simulations. 

4.1 Results for relative consistency 
4.1.1 Rank one initial density matrix 

Fig. (0) shows the probability tree and minimum consistency statistics for a simu- 
lation with e = 0.03 ~ Q.2/\fd and a rank one initial density matrix. As expected 
there is one almost immediate trivial projection. Three more projections occur and 
then no more. From fig. (|l|) the probability of a projection with five histories and 




(a) probability tree (b) consistency statistics 

Figure 3: relative consistency, e = 0.03 ~ 0.2/ y/d 

e = 0.03 is less than 1% so this result is as expected. The simulation was run for 
longer than is shown in the figure (until t = 100) but no further projections occurred. 

Fig. (Q) shows the probability tree and projection times for a simulation with 
e = 0.15 « 1/y/d. A gain there is the initial trivial projection but no other trivial pro- 
jections occur. Projections then occurred at roughly equal equal time intervals until 
there were fifteen histories. The time between projections then rapidly increased. 
This is in accord with fig. ([!]) as the probability for a projection with fifteen histories 
and e = 0.15 is around 5%. Projections after this time only occur for large deviation 
away from the mean and therefore occur extremely erratically. These later projec- 
tions are extremely unlikely to vary smoothly for a range of e. The simulation was 
run until t « 100 and no further projections occurred. This simulation has produced 
an interesting set of histories with a complicated branching structure. 

The next pair of figures fig. (||) shows the results of a simulation with all the 
parameters unchanged except for e which is now 0.16. The qualitative description is 
the same and the first eight or so projections are similar. After that however the two 
sets of histories are very different. This is the problem with the algorithm applied to 
this model: interesting sets of histories are produced, but they change dramatically 
for small changes in e. 

From fig. (|I|) choosing e = 0.3 ~ 2y/d looks large enough so that projections will 
always be made before the background level is reached. The theoretical analysis 




(a) probability tree (b) projection times 



Figure 4: relative consistency, e = 0.15 ~ 1/y/d 

also suggests that for such a large value of e some repeated projections will occur. 
Indeed fig. (|6]) demonstrates that nine repeated projections occurred each giving a 
history of probability S. 

An interesting alternative is to choose e as a percentile from fig. ([]]), that is e(k) pa 
e(e) \ogk where k is the number of histories. Fig. (|7|) demonstrates the consistency 
statistics for a run with e(k) chosen at the 50% level. All of the probabilities except 
for the initial projection were non-trivial. Rather than the projections being made 
in regimes where the DHP fluctuates about its mean value most of the projections 
have been made at times when the DHP is monotonically decreasing, so that the 
histories are much more likely to vary continuously with e. Two other advantages of 
choosing e this way are that larger sets of histories are produced, and if an algorithm 
is designed to produce a set of histories of a certain size choosing e in this way will 
produce a more consistent set than choosing e to be constant. However, though the 
results are more stable (when the percentile is changed) than for constant e, results 
from simulations show that they still change too much to single out a definite set of 
histories. 

By looking at the consistency statistics the problem is easy to understand. Since 
a projection is made at the earliest possible time generically once it has been made 
the DHP jumps up as the most consistent projection has occurred. The consistency 



0.0000 

(a) probability tree 



projection number 

(b) projection times 



Figure 5: relative consistency, e = 0.16 ~ 1/Vd 
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Figure 6: relative consistency, e = 0.3 ~ 2/\fd 

level then falls. While it is decreasing monotonically any change in e will produce a 
continuous change in the time of the next projection. However, if e is too far below 
its mean level the projection times will vary discontinuously, and all the projections 
afterwards will generically be completely different. Since the mean level of the DHP 
depends on the number of histories strongly either e must be chosen sufficiently large 
so as to be above this or it must be chosen so as to increase with the number of 
histories. This is demonstrated by the first few projections as shown in fig. (§) — a 
close-up of fig. (0b) would also show this. 
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Figure 7: relative consistency, e chosen at 50% 



4.1.2 Other initial conditions 




(a) probability tree (b) consistency statistics 

Figure 8: relative consistency, 5 = 0.02, e = 0.05 

Simulations with a rank 2 initial state, and all other parameters remaining the 
same, produce the same results except that each of the initial projections is repeated 
producing two trivial histories with probability 5. We can choose 5 according to 
eq. ( p2"D to try to prevent these projections, that is choose 5 = 0(e 2 ). Since many 
of the histories we expect to generate will have probabilities smaller than this it is 
sensible to use a relative non-triviality criterion^. The analysis leading to eq. (|22|) is 

4 Using a relative non-triviality criterion earlier does not qualitatively change the results except 
that the initial trivial history would also have been extended — the results would have been 
qualitatively the same as for rank two initial reduced density matrices. 



only accurate to first order in t, therefor eq. (p^ ) is only valid when 8 is sufficiently 
small. If 8 is too large the consistency level of a reprojection will start decreasing 
and when a reprojection eventually becomes non-trivial it will be consistent. For 
the example discussed there were no values of 8 with e = 0.15 that prevent an initial 
trivial reprojection. Fig. ([8]) demonstrates the start of a simulation with e = 0.05 and 
8 = 0.02 = 8e 2 . The graph of the consistency statistics shows the initial projection at 
t = and then that there are no non-trivial extensions until t ~ 0.08 by which time 
the projection is not consistent. A non-trivial projection is made at t ~ 0.32 and the 
algorithm then continues as before, with the trivial projection avoided. Because the 
projection has not occurred with probability 8 there is a range of values for 8 that 
do not affect the resulting histories — they are independent of 8. The first three 
projection times will obviously vary continuously for a small range of e. The other 
projections that occurred in this simulation all occurred at much more separated 
times in a regime where the consistency level was not decreasing monotonically. If e 
is chosen according to the percentile distribution in fig. (JT|) — that is e(k) ~ e(e) log k 
where k is the number of histories — e will be small enough initially to prevent any 
trivial projections (with an initial density matrix of rank greater than one) and will 
allow a full set of histories to be built up at later times. 

Simulations with a full rank initial density matrix and e = 0.15 result in a trivial 
repeated projection for each initial history. If e is smaller (< 0.07) no trivial repro- 
jections occur. If two initial projections are made uncorrelated with the Schmidt 
states then further (Schmidt) projections at t = will not generically be trivial or 
consistent. In both cases, after the initial projections and possible reprojections the 
qualitative behaviour is the same as the rank one case. 

4.2 Results for absolute consistency 

To produce interesting sets of histories from an absolute consistency criterion three 
effects need to be balanced against each other. If e is too small the most likely 
projections will be those that produce very small probability histories. If 8/e is too 
small the likelihood of repeated projections (hence trivial histories) will be high. 
If 8 is too large the non-triviality criterion will dominate the algorithm and only 
probability 8 (trivial) histories will be produced. Only an absolute parameterised 
non-triviality criterion is considered since a relative criterion will clearly produce 



.01 

Figure 9: absolute consistency, 5 = 0.01, e = 0.1 

almost nothing except infinitesimal histories. The following results demonstrate 
these effects. 

4.2.1 Rank one initial density matrix 

For example if 5 < e and the initial reduced density matrix has rank one the algo- 
rithm will generically produce \1/S\ trivial histories. This is a particularly simple 
case of the analysis that suggests that if 5/e < 0(1) repeated projections are prob- 
able. Fig. (|9]) shows an example of this from a computer simulation with e = .1 
and 5 = 0.01. Only the first ten projections are shown. This behaviour remains the 
same in the limit as e — > 0, 5 — > 0, 5 < e. 




probability tree 
Figure 10: absolute consistency, 5 = 0.01, e = 0.02 

As the ratio 5/e increases and becomes 0(1) the nature of the set of histories 
changes. Occasionally when a reprojection becomes non-trivial it will no longer be 



consistent and a reprojection will not occur. A significant time may elapse before the 
next projection is made which will result in a non-trivial projection, which will then 
be followed by more trivial repeated projections. This is demonstrated in fig. (fCDD 
where e = 0.02 and 5 = 0.01. Though this is an interesting set of histories this range 
of parameter values does not give a theory with predictive power since simulations 
show that the results vary enormously for small changes in e and S. 




(a) probability tree (b) consistency statistics 



Figure 11: absolute consistency, 5 = 0.01, e = 0.001 

As S/e increases past 1 the number of histories made with probability 5 decreases 
to just the initial projection. Fig. flTT| ) shows that this range of parameter values 
produces interesting histories but the projections are occuring at times when the 
consistency level is fluctuating randomly about the mean and so will be unstable to 
small changes in e. 

4.2.2 Other initial conditions 

Choosing larger rank initial reduced density matrices or initial projections does 
not qualitatively change the analysis. The only difference is that for 5 > e and e 
sufficiently small no trivial projections will be made. 



5 Conclusions 



The algorithm produces sets of histories with a complicated branching structure 
and with many non-trivial projections for a range of parameter values. Algorithms 
using the DHC produce results that are essentially the same for a wide range of S 
(the non-triviality parameter) including the limit 5 —>■ 0. However, the algorithm 
does not make useful predictions when applied to this model since the results vary 
erratically with e and there is no special choice of e singled out. Choosing e as a 
function of the number of histories according to fig. (jl]) produces the least unstable 
sets of histories and the largest sets of non-trivial histories, but even in this case the 
algorithm does not single out a definite set. 

The algorithm is less effective when used with an absolute consistency criterion: 
in this case the predictions of the algorithm also vary erratically with 5 and the 
resulting sets of histories include fewer non-trivial histories. 

The results of the simulations agree well with the theoretical analysis of sec- 
tion (Q) and demonstrate features of the algorithm that will also apply to other 
models — such as the analysis of repeated projections. They also demonstrate some 
of the difficulties that an algorithm must overcome. These problems can be related 
to the discussion of recoherence in refs. ||. The algorithm will only produce stable 
results (with respect to e) if the projections occur when the off-diagonal terms of 
the decoherence matrix are monotonically decreasing. This behaviour is only likely 
in a system like this for times small compared to the recurrence time of the system 
and when the number of histories is small compared to the size of the environment 
Hilbert space. The results of the model do show stability for the first few projections 
and if much larger spaces were used this behaviour would be expected for a larger 
number of histories. In particular as the size of the environment goes to infinity it 
is plausible that the algorithm applied to this model will produce a large, stable, 
non-trivial set of histories. 

In retrospect it was ambitious to hope that an algorithm applied to this model 
would produce large sets of stable histories. A random model like this, with such a 



small environment, will generically only decohere a few histories (see also [14, 15 



The Gaussian Unitary Ensemble 



In the GUE the matrix elements are chosen according to the distribution 

2 n/2 f TrM 2 )l 

P( A ) = r/n-M/s-^ ex P 



[(2tt) 1 /V]« 2 ^ [ Aa 2 

2"/ 2 ^ f .4, ,.4/, 
[(2tt) 1 /2 (T ] 



n>j,fe>l 



II expj-^g) [] expj-^g) J] exp^ - ■■>'- 



LV Z7i 4 a J n>j>=l I ^ a J n>k>j>l { ZIJ ) n>k>j>l { ZIJ ) 

where Aj k = Xj k + iYj k , Xj k = X k j and Yj k = —Y k j. Therfore all the elements are 
independently, normally distributed, the diagonal with variance la and the real and 
imaginary off-diagonal with variance a. Some expectations for a normal variable 
with variance a are 



E[\X\ n ] 



2 „/ 2cT „ r (£±i) 

and in particular E[\X\] = \/(2/n)a, E[X 2 ] = a 2 and E[X 4 ] = 3a 2 . Since X tj is 
independent of X k i unless i — k and j = I, or i = I and j — k 

E[XijX k i] = a 2 (5 u 8 jk + S ik 5ji), 
E[YyYiw] = a 2 (5 ik 5ji- 8 a 5 jk ). 

Therefore, for the elements of A 

E[Aj] = o, 
E[^ij>l fc i] = 2a 2 5 u 5 jk , 
E[Ay A k [A mn A op ] = 4(7 (SuSj k 5 mp 5 no + 5i n 5j m 5 kp 5i + 5i P 5j 5 kn 5i m ). 

Applying these results to vectors and projection operators, 

E[nUm] = 0, (36) 
EflnMrnl 2 ] = 2 ( x 2 |n| 2 |m| 2 , 
E[n t 4P4m] = 2Ax 2 n t m, 



EflnUPAml 2 ] = 4a 4 [d 2 |n t m| 2 + (n t Pn)(m t Pm) + d|n| 2 |m| 2 ], 

where d is the rank of P. For the real part or the imaginary part just take half of 
the above since \z\ 2 = [Re(z)] 2 + [Im(2:)] 2 . 
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